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Inflationary observables, like the power spectrum, computed at one- and higher-order loop level seem 
to be plagued by large infra-red corrections. In this short note, we point out that these large infra- 
, red corrections appear only in quantities which are not directly observable. This is in agreement 

^ |. with general expectations concerning infra-red effects. 

P , PACS numbers: 98.80.Cq CERN-PH-TH/2007-214 

^ ' Cosmological inflation [Ij has become the dominant paradigm within which one can attempt to understand the initial 
conditions for Cosmic Microwave Background (CMB) anisotropics and structure formation. In the inflationary 
' picture, the primordial cosmological perturbations are created from quantum fluctuations which are "redshifted" 
^ , out of the horizon during an early period of accelerated expansion. Once outside the horizon, they remain "frozen" 
' until the horizon grows during a later matter- or radiation-dominated era. After falling back inside the horizon they 
\ are communicated to the primordial plasma and hence are directly observable as temperature anisotropics in the 
. CMB. These anisotropics were first detected by the Cosmic Background Explorer (COBE) satellite [1,11,13, ^^'^ 
more recently they have been mapped with spectacular accuracy by the Wilkinson Microwave Anisotropy Probe 
(WMAP) 11 . 

J — These CMB observations show that the cosmological perturbations are very small, of order 10^^ compared to 

, the homogeneous background. Therefore, one might think that first-order perturbation theory will be adequate for 
' all comparison with observations. However, that may not be the case; the Planck satellite [6] and its successors 
. !^ \ may be sensitive to non-Gaussianity (NG) in the cosmological perturbations. Such non-Gaussianities are sourced 
• by self-interactions in the early universe, and become visible at the level of second- or higher-order perturbation 
5—1 ■ theory This possibility is of considerable interest in its own right and is presently being vigorously explored. 
However, there are other equally compelling reasons to go beyond linear theory. Self-interactions of any scalar field 
during the inflationary stage, and more interestingly of the comoving curvature perturbation C, do not only imply 
non-Gaussian statistics: they also give rise to corrections in all correlators, including the observationally interesting 
cases of the two- and the three-point correlation functions. Such corrections are associated with so-called "loops," 
in which virtual particles with arbitrary momentum are emitted and re-absorbed by the fields which participate in 
the correlation function. 

Loop corrections may lead to significant effects [1, S [13: [HI 113 which scale like (powers of) the number of e- folds 
between horizon exit of the mode k under consideration and the end of inflation. Furthermore, it has been known for 
a long time that there exists a breakdown in the perturbative expansion due to infra-red (IR) divergences [isl . . 
These classical one-loop corrections depend on the total number of e-folds of inflation since they scale like powers 
of In a, where a is the scale factor of the universe [l^ [H, [13, [H, [H, [lO] . 

To illustrate the appearance of IR divergences, we can adopt a simple toy model made of a scalar field cj) with 
cubic interaction ^(j>^ This field is not necessarily the scalar field driving infiation. By using the in-in formalism 
suitable for non-equilibrium field theory, we can compute the one-loop correction to the power spectrum of the 
scalar field perturbations in a pure de Sitter epoch characterized by a constant Hubble rate H. The computation 
is performed in Appendix A, to which the reader is referred for all technical details. One finds 
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where r is the conformal time and L^^ is the infra-red comoving momentum cut-ofF. One is now faced with the 
task of determining what value should be assigned to L~^. 

One option is to set L to be some comoving length scale which left the horizon many e-folds before the observable 
universe. In particular, the smallest possible value of is a^TJ, where is the value of the scale factor at the 
beginning of inflation. Since the wavelength goes out of the horizon when the scale factor equals = k/H, we 
see that ln{kL) = ln(afe/ai), which is proportional to the total number of e-folds from the beginning of inflation to 
the time when the mode k exits the horizon. This is the large cumulative IR correction previously alluded to: it gets 
larger the longer inflation lasts. Furthermore, the In^(fcr) accounts for the the number of e-folds between horizon 
exit of the mode k and, eventually, the end of inflation. More generally, at the generic n-th order of perturbation 
theory, the power spectrum is expected to get corrections growing like ln"(/cL). 

Faced with these potentially large logarithmic corrections, one may take the road of trying to resum them, for 
example by using the standard method of the Wilsonian Renormalization Group (RG). This technique, however, 
may not be that efficient in practice. Indeed, the full set of exact RG equations need to be solved. Approximations 
inevitably lead to disregarding a set of diagrams which are not at all subleading in the IR. We refer the reader to 
the Appendix for further discussions about this point. 

More relevantly, before attacking the problem of the large IR corrections, we should first ask ourselves whether 
they are really present in any physical observable which can be measured. Indeed, one should be suspicious that 
large effects from infra-red modes are somehow unphysical, because on general grounds we expect the structure 
within our observable patch of the universe to depend only on its local properties. 

What we would like to point out in this short note is that the large IR corrections appear only in quantities 
which are not directly observable. Although their appearance cannot be avoided, they do not themselves have any 
particular interest unless one wishes to ask more general questions, for instance, what is the probability to find an 
inflaton field homogeneous enough to lead to the correct temperature anisotropy in the CMB within patches like 
ours? Their presence is an indication that we have computed the answer to a question which we can never observe. 
When large IR corrections appear, we should instead learn how to change the questions we ask in order to obtain 
predictions which are genuinely local and apply within our observable patch. Some of these considerations already 
appeared in Refs. p^ . [20I [23j and in the present work we further elaborate them. 

Let us consider the primordial curvature perturbation computed using the 6N formalism [U [25l [26l . [27| , where 
N is the number of e-folds from an initial flat hypersurface to a final uniform-density hypersurface. It is convenient 
to pick the initial time to be shortly after Hubble-exit of the scale over which one is smoothing the relevant physical 
quantities. The final slice has to be located before the smoothing scale re-enters the horizon in such a way that 
the separate universe approach may be employed [1^ [2^, [lO, HH . The homogenous Fricdmann- Robertson- Walker 
(FRW) equations can be used to calculate N on large (super-Hubble) scales. 

Assuming for simplicity the presence of a single scalar field 0, the inflaton, the primordial curvature perturbation 
C is determined by the value of the scalar field on the initial slice 

C(x) = 5iV(</)(x)) = 7V(^(x)) - N(^) = N'ScPi^) + iiV"j02(x) + . . . , (2) 

where 

and so on. (3) 

The connected two- and three-point functions of the scalar field are defined by 

(0ki0k.) = F0(fc)(27r)3<53(ki+k2), 
(0ki</'k.0k3) = B4fci,A:2,fc3)(2^)^'5'(ki+k2+k2), (4) 

where ^^(fc) is related to the power spectrum 'P^{k) by 'Pip{k) = {k'^ P^{k) / 2tt'^) . The presence of a nonvanishing 
bispectrum B^(ki,k2,kj,) accounts (at this order) for the possibility that the scalar fiuctuations are not Gaussian 
at horizon exit and, for simplicity, we neglect the trispectrum and higher-order connected functions. 

We will restrict ourselves from now on to the power spectrum of the primordial curvature perturbation, C,, 
which is a physical observable within our local patch. In order to obtain predictions concerning its properties, all 
computations should be done within a comoving box whose present size I is not much larger than the present horizon 
in order that the answer depend only on local properties of the universe at our position. In the language of 
Ref . [23| , this may be called a minimal box. We shall also be obliged to consider a superlarge box, whose comoving 
size L left the horizon right at the beginning of the inflationary stage. What we would like to stress in the following 
is that the appearance of large logarithmic IR divergences is a consequence of adopting the superlarge box instead 
of the minimal one [2^ . 
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We have said that in order to make predictions for we must compute within a box not much larger than . 
However, although we are free in principle to carry out the computation within a box of any convenient size, in 
practice we may be restricted by our inability to compute the relevant correlators. For example, in chaotic inflation 
beginning at energy densities close to the Planck scale, the original inflationary patch will be exponentially larger 
than the horizon scale when modes corresponding to the CMB left the Hubble radius. In general, our only choice 
at present is to compute correlation functions within the entire inflationary region, which constitutes a superlarge 
box. However this computation need not have anything to do with the anisotropy in the CMB temperature: the 
correlation functions within the superlarge box are essentially averages of the correlation functions computed within 
a horizon-sized box, where the average is taken over all ways one can fit the small box within the superlarge one. This 
correlation function will only successfully predict the CMB anistropy if conditions in our universe yield correlators 
which are very close to those computed by averaging over the superlarge box. 

A special simplification occurs if the universe contains only a single field. In this case, one knows the conditions 
leading to the end of inflation and (providing the background field is sufficiently homogeneous) one therefore always 
has the choice to compute in a horizon-sized box. However, where isocurvature fields are present it will generally 
be true that these fields assume different values in different regions of the superlarge box. (This issue is discussed 
in more detail below.) In that case, one does not know a priori the conditions leading to the end of inflation in 
each small box, and the only alternative in analytic calculations is to compute in the superlarge box instead. 

What is the relationship between the correlation functions in differently sized boxes? As one varies the small box 
within the superlarge one, there will be a slow variation in background quantities such as the Hubble parameter 
H, the slow-roll parameters e, 77, ^ and so on, together with all other zero-momentum quantities in the theory^. 
The expectation value in the superlarge box is found by taking averages over these background quantities. If the 
background quantities exhibit large variations within the superlarge box then these averages will develop significant 
contributions in the infra-red, associated with wavenumbers between and the inverse size of the small box. 
On the other hand, we could equally well have obtained the expectation value in the superlarge box by computing 
the ensemble average from first principles, and using the ergodic theorem to connect this with the spatial average. 
These two methods of computation must agree. Therefore, we see that large IR divergences are merely associated 
with large fluctuations in the background fields on scales very much larger than the presently observable region. 
These fluctuations may be highly non-gaussian, but because they are restricted to scales enormously larger than 
any which are accessible to experiment, they are of no observational interest. The presence of such structure on 
ultra-large scales was anticipated long ago by Salopek & Bond [29| . 

To strengthen the conclusion that the average {Vi() coincides with the power spectrum 'P(^ computed within the 
superlarge box, we follow Ref. [23| and consider a box of generic size M <C i placed within the superlarge box. We 
want to show that a quantity such as {Pm(), and consequently the average of the power spectrum {Vmc), does not 
depend upon the scale M and is therefore a reasonable candidate for the expectation value within a box of size L. 

The curvature perturbation within the box of size M is given by the 5N formula, which yields 

C(x) = n'.Mm (x) + \N'j4>lt (x) + • • • , (5) 

where now 

and so on. (6) 

We wish to relate this to the same computation performed within the superlarge box, which entails replacing all 
quantities evaluated within M by their equivalents evaluated within L. This can be accomplished by using the 
separate universe picture to account for variation in the background field, which gives 

5(t)M = (0- 0m) + '^'/'(x), 

n'm = n'^ + N'i (0 - 0,,) -f iiv;" (0 - 0,,) % . . . . (7) 

For the sake of simplicity, let us focus on the contributions to the two-point correlator coming from the bispectrum 
(if the scalar field pertubations are Gaussian at horizon exit, the proof is contained in Ref. |23i]). From Ref. (32] we 



(x) = 0(x) - 



dN 



This sounds contradictory, but such quantities are to be thought of as "zero-momentum" only within the small box. Within the large 
box they are a superposition of modes with wavenumbers in the far infra-red. When moving to some larger box, one redefines the 
spatial average of such a quantity as its new "zero-momentum" value, and shifts some of the infra-red region into the perturbations. 
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may directly read off the contributions to Pmc from the bispectrum. They start at one-loop and up to two-loops 
they read 

Pmc""" ^ I NmN'mBm4> (k, q, \k, - q|) , 

PmI""" 3 -77^ I d^qi d\2 nIiNZBm^ (|ki + qi| , |q2 - Qil , \\^2 " q2|) PmM)^ (») 

where only momenta larger than AI^^ are to be considered because the box introduces periodic boundary conditions, 
and BmiP is the bispectrum of the inflaton within the M-sized box under consideration. These contributions are 
accompanied by terms coming from the power spectrum, as described in Ref. ^32J . which we are ignoring here. 

In any M-box, the tree-level contribution to Pmq is (-^m)^-Pm0- Within the L-sized box we can use the separate 
universe picture to write 

(N'm^Pmc^ = [N'^fPL^ + {{N'lf + N'^N'l' + (2e + r,){N'^f + AV2eN'^N'l] {^,, - ^Pl^ + ■■■ , (9) 

where a term linear in {(pj^j — cj)) has been omitted and e — (mp[/47r)(iJ'(0)/i?((/)))^, 77 — {m^^/ At:)H" {(j)) / H{(f)) 
are the usual slow-roll parameters. The term proportional to (-/V£)^ -I- N'j^N'l' cancels that part of the running with 
M which is sourced by the omitted power spectrum terms in Eq. ([8]), as described in Ref. [2^. Of the remaining 
terms, (2e + rj){N'^Y is of the same functional form (in derivatives of N) as the leading term {N'j^)'^Pl^, but is 
suppressed by extra slow- roll parameters. This term will cancel contributions from subleading terms in which 
are not visible when we compute only to leading order. Therefore this term can be ignored. The interesting term is 
the one proportional to 4:^/2€N'j^N'l, which is of the same functional form in derivatives of N as the one-loop term 
in Eq. ^. 

One must average Eq. ^ over all possible positions of the AI -sized box within the large box. S mce — J ^s 
simply 5(j}{n) (within the large box) smoothed with a top-hat window function, its average is given by 

((0M-'^)'>= j^^ \v<,{k). (10) 

Now consider how the one-loop term in Eq. ([8]) runs with AI . Within the L-sized box, one can write (in analogy 
with Eq. ^) 

Pit"" ^j:^ I {N'^N'^Bl^ + S(0M - 0)' } + • • • , (11) 

where a linear term has again been discarded, and B is an abbreviation for the combination 

B ^ [^N'lN'l' + In'.K'^ Bl^ + {N'^N'l' + (TVi')'} + \n'lN'IBI^. (12) 

In this equation, B'j^^ and _B^^ are, respectively, the first and second derivatives of the inflaton bispectrum within 
the i-sized box with respect to the background field, and all the bispectra are evaluated with arguments fc, q and 
|ki — q|. The leading term in Eq. pT|) is proportional to N'j^N'IBl,^- The running of this term comes entirely from 
differentiating the integral with respect to its limit. Using the method of Boubekeur & Lyth [s^ to estimate the 
integral, and the squeezed bispectrum computed by Allen et al. [s^l (see also Ref. [1^), it follows that this running 
exactly cancels the running from the term proportional to N'j^N'l in Eq. 

The remaining terms in Eq. ()11|) . which are described by B, acquire running from both the lower limit of the 
integral, and the background expectation value {{4>m ~ They are proportional to three copies of P^ and 

therefore are formally comparable to two-loop corrections. In a full calculation, one should include all two-loop 
terms and show that the running of the B terms correctly cancels that of the two-loop contributions. In this paper 
we focus on the two-loop term written in Eq. ([5]) and show that its running cancels with the term of the same 
functional form (in derivatives of N) in B. The running imparted to this term by {{4>f^j — </>)^) is given by 

^(Pl7c'°°'^) 2 ^^iV>;>i^(A/-l) j^^ ^ d?q Bl^ (fc, q, |ki - q|) . (13) 

This can be thought of as a renomalization of the classical value of the scalar field as seen inside the small box. 
On the other hand, the Af-dependence of P^.j^°°^ comes by differentiating the integral with respect to its limit 
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where we have assumed that the relevant physical momenta are much larger than AI^^. The contribution (|14p cancels 
two thirds of the running in Eq. ^3]) . The remaining part is associated with "dressing" of the 6N coefficients, as 
described in Ref. [s^l- It is cancelled by a bispectrum term coming from the N'^N'^' term in the 5N expansion, 
leaving zero running overall. A similar computation may be performed including the trispectrum. 

The computation above confirms that {Vmc) does not depend on the size of the box M, as expected on general 
grounds, and that it coincides with the power spectrum computed in the superlarge box: large IR divergences 
inevitably appear because what we are computing is in fact (7^£^), that is the power spectrum on the superlarge 
box. 

One can choose a small box to minimize the two-loop contribution, but this makes the one-loop contribution 
sensitive to the largest scale L. For instance, in slow-roll models of inflation = 0{e^^^)P?. It can be shown that 



which becomes significant for M <^ L. Alternatively, one can try to reduce the one-loop contribution going to the 
superlarge box, paying the price of a large two-loop contribution. The point though is that the large IR divergences 
show up in the average power spectrum and, as we have advocated before, this is not a quantity of physical interest. 
It only provides the level of uncertainty inherent in the theoretical predictions. 

If one insists in adopting a superlarge box, it should be kept in mind therefore that the location of our box may 
be untypical and one should quantify how likely it is that the correlators averaged over the superlarge box coincide 
with the correlators in our observable universe. To deal with this problem, one can try to use the approach of 
stochatic inflation. This embodies the idea that the IR part of the scalar field may be considered as a classical 
space-dependent stochastic field satisfying a local Langevin-like equation [36i |. The stochastic noise terms arise 
from the quantum fluctuations which become classical at horizon crossing and then contribute to the background. 
One can derive a Fokker-Planck equation describing how the probability of scalar field values at a given spatial 
point evolves with time. One should split the computations into two steps. First, the distribution of the values 
of the scalar field inside the superlarge box should be estimated. This will allow one to compute the probability 
P(0gQ I igo) that the starting initial condition for the background field is given by = (j)QQ in our local observable 
patch at time tgOi j-e. when there remain 60 or so e-folds till the end of inflation in that patch. The distribution 
of the scalar fleld in the superlage box, due to the IR divergences, will be far from being Gaussian. However, we 
do expect that in our local patch initial conditions are relatively uniform. From the point of view of our local 
universe the IR properties of the superlarge box only enter to provide the probability P{4>qq \ ieo) that the value 
(/)go is achieved in our patch. At later times the field value at a given point has a distribution given by 



Gaussian distributed as large IR corrections are present in the superlarge box [29|, |37| . However, if we impose the 
observable constraint in order to select those initial parameters which lead to CMB anisotropies consistent with 
observations, a severe selection is imposed on initial conditions and we expect a probability close to Gaussian. In the 
stochastic approach, therefore, one trades the uncertainty in the prediction coming from the large IR contributions 
in the superlarge box for the uncertainty inherent in having a probability distribution for the background quantities. 
It would be interesting to determine such a probability distribution to evaluate the uncertainty in the theoretical 
predictions for our local universe. 
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APPENDIX A: THE IN~IN FORMALISM FOR NON-EQUILIBRIUM QUANTUM FIELD THEORY 

In this Section we briefly present some of the basic features of the non-equilibrium quantum fleld theory based on 
the in-in, also dubbed Schwinger-Keldysh, formulation '21]. The interested reader is referred to the excellent review 
by Chou et al. ^] for a more exhaustive discussion. 
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FIG. 1: The appropriate time contour C. 



Since we need the temporal evolution of quantum correlators with definite initial conditions and not simply the 
transition amplitude of particle reactions, the ordinary equilibrium quantum field theory at finite temperature is not 
the appropriate tool. The most appropriate extension of the field theory to deal with nonequilibrium phenomena 
amounts to generalizing the time contour of integration to a closed-time path. More precisely, the time integration 
contour is deformed to run from — cx) to +00 and back to — cx3. 

The CTP formalism (often dubbed as in-in formalism) is a powerful Green's function formulation for describing 
non-equilibrium phenomena in field theory. It allows to describe phase-transition phenomena and to obtain a 
self-consistent set of quantum Boltzmann equations. The formalism yields various quantum averages of operators 
evaluated in the in-state without specifying the out-state. On the contrary, the ordinary quantum field theory 
(often dubbed as in-out formalism) yields quantum averages of the operators evaluated with an in-state at one end 
and an out-state at the other. 

For example, because of the time contour deformation, the partition function in the in-in formalism for a real 
scalar field is defined to be 



Z[J] = Tr 
= Tr 



X T_ 



T l^exp 
T+ ^exp 
exp 



i / J(j) 
Jc 



(Al) 



where C in the integral denotes that the time integration contour runs from minus infinity to plus infinity and then 
back to minus infinity again, see Fig. [TJ The symbol p represents the initial density matrix and the fields are in 
the Heisenberg picture and defined on this closed time contour. Sometimes it is more usful to work with other 
field variables, (f>c = 1/2(0+ + </)_) and (pA ~ {4>+ ^ </'-)■ In such a case one has to properly redefine the sources 
as Jc = 1/2(J+ -I- J_) and Ja = (•/+ — J-)- To identify the physical degrees of freedom, the normalization of the 
generating functional Z [J a, Jc]\j^=o = 1 has to be imposed. 

We must now identify field variables with arguments on the positive or negative directional branches of the time 
path. This doubling of field variables leads to six different real-time propagators on the contour [l^j. These six 
propagators are not independent, but using all of them simplifies the notation. For a generic bosonic neutral scalar 
field 4> they are defined as 



G-+ix,y) 
G+- {x,y) 
G++ {x,y) 
G~-{x,y) 
G«(a;,y) 
G^ix,y) 



i{(t){x)(t){y)), 

i{(f){y)(l){x)) , 

G-+{x,y)9{x,y) 

G+-{x,y)e(y,x) 

G++ix,y)-G+- 

G++(x,y)-G-^ 



^G+-{x,y)eiy,x), 
^G-+{x,y)9{x,y), 
ix,y) = {G-+ix,y)-G^ 
ix,y) = {G+-{x,y)-G- 



'{x,y)) 0{xo - yo), 
'{x,y)) 0{yo - xo), 



(A2) 



where the last two Green functions are the retarded and advanced Green functions respectively and 9{x,y) = 
9{xo ^ yo) is the step function. When computing a loop diagram, one has to assign to the interaction points a plus 
or a minus sign in all possible manners and sum all the possible diagrams, taking into account that vertices which 
a minus sign has been assigned to must be multiplied by —1. 

In the basis of the fields 4>c and c/ja, one may define the Green functions 



{M^)Uy)) = G'^{x,y) = --{G-+{x,y) + G+-{x,y)), 

{M^)4>A{y)) - -*G^(x,y), (A3) 
while the Green function {4>A{x)(j3A{y)) vanishes identically because of the identity G^ + G^ = + G ^ . 
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FIG. 2: The cubic vertices with the proper signs + and 




FIG. 3: The diagram contributing to the one-loop two-point correlation function. To each vertex the signs -I- and — should 
be attached to each vertex according to the correlator one wishes to compute. 



For equilibrium phenomena, the brackets (• • • ) imply a thermodynamic average over all the possible states of 
the system. While for homogeneous systems in equilibrium, the Green functions depend only upon the difference 
of their arguments {x,y) = {x — y) and there is no dependence upon {x -\- y), for systems out of equilibrium, the 
definition (|A2[) has a different meaning. The concept of thermodynamic averaging is now ill-defined. Instead, the 
bracket means the need to average over all the available states of the system for the non-equilibrium distributions. 
Furthermore, the arguments of the Green functions (x, y) are not usually given as the difference [x—y). For example, 
non-equilibrium could be caused by transients which make the Green functions depend upon (a;o,?;o) rather than 
{xo - yo). 

The Lagrangian we consider is the one for a massless (at the tree level) scalar field with a cubic self-interaction. 
This field is not necessarily the scalar field driving inflation. 



C[cb+] - = i^-g"Pid,cb+dp<P+ - d.c^^dpcj^^) ~ !^^ {4>l - ) j ■ (A4) 

We are using the conformal metric g^p — o?{t) diag(I, —1, —1, —1), where a(T) = — l/iJr (r < is the conformal 
time) is the scale factor during the de Sitter stage characterized by a Hubble rate H. The interaction term is cubic 
in the scalar field, see Fig. ^ and /i has therefore the dimensions of a mass. 

Since we will be interested in the following in the infra-red divergences, we provide here the expressions for 
the Fourier transformed of the tree level expressions g'^ and g^ of the two-point correlation functions and 
respectively in momentum space for wavelengths larger than the horizon 



^2) 0^ fff (ri, T2) c ^0{t, - T2) [k^T^ - r|)] , (fcn, fcr2 « 1). (A5) 

As an illustrative example of the infra-red divergences, we compute the two-point correlation function G'£(t, r'). 
The relevant diagram is given in Fig. [3] which represents the sum of a set of diagrams whose relative sign is dictated 
by the appropriate sign of the cubic vertex. 

A straightforward calculation which makes use of the expressions (|A5P , shows that the infra-red contribution to 
G|(t, t') at one-loop is given by 



FIG. 4: The diagrams contributing to the two-loop two-point correlation function. To each vertex the signs + and — should 
be attached to each vertex according to the correlator one wishes to compute. 



+ J C?na^(Tl) J dT2a\T2)g^{T\Ti)gl{T,T2) J J^9i+p(.ri,T2)g^{Tl,T2) 



9 (2^)2^2 yj 'J ' {T1T2Y 'J '^'^ J {tiT2 



{T1T2Y / Wl-1 P 



(A6) 



where is the comoving infra-red cut-off which we may set to be equal to Oii/, being the value of the scale 
factor at the beginning of inflation (the upper limit of integration over the momenta p is dictated by the fact that 
the expression is valid for momenta p <k). Notice that the power spectrum 7^^ of the perturbations of the scalar 
field (j> is directly related to the correlation function G'^{t,t') computed at equal times 

V4k) ^ ^GI{t,t). (A7) 
A similar computation leads to the one loop correction to the retarded Green function in the infra-red 



Gf-'°-(r,r') ^ I 'liyH^'^ HkL)lAr/r'). (A8) 
At two-loop, one can easily calculate the IR contribution of the first three graphs in Fig. [4] obtaining 



gfl,2-ioop, ,^ _ g^(r, r') In^ffcL) In^ t/t'). (A9) 

It is not difficult to show that the same kind of topology of diagrams appearing at higher loops leads to the 
resummation 



Gf-(r,r') ^ 5f (T,r')exp (^^^^-^ In(fcL) ln^(r/r')) . (AlO) 

However, diagrams like the forth one in Fig. (U) are not automatically included in this resummation. This is 
unfortunate, because they are not subleading. This is reminiscent of what happens in a At/)"* theory where, performing 
a Hartree approximation, one is able to resum the all class of "daisy" diagrams; however, unless an 0{N) symmetry 
is invoked for large N , one has no guarantee that the omitted diagrams are unimportant. 



d/dA =1/2 



.ci 



FIG. 5: The RG equation for the two-point correlation function. 

In general, one can hope to resum the large leading logarithms using the standard method of the Wilsonian 
Renormalization Group (RG) techniques developed in quantum field theory and in statistical physics [ssj . These 
techniques have been recently applied to resum cosmological perturbations of a coUisionless fluid ^9] . The starting 
point of this formulation of the Renormalization Group is a modification of the generic two-point correlation function 
as follows 

Gp(ri,T2) -> G^{n,T2) = Gp{n,T2)9{p,A), (All) 

where 9{p, A) is a high-pass filtering function, which equals unity if p 3> A and zero if p ^ A. The modified 
generating functional Z\[J] can be obtained from the expression (jAip and describes a fictitious universe in which 
all fluctuations with momenta smaller than A are damped out. In the limit A —> all fluctuations are included 
and we recover the physical situation. Decreasing the cut off, the linear and the nonlinear effects of smaller and 
smaller fluctuations are gradually taken into account. This process is described by the RG equation obtained by 
taking the A derivative of Za 



where the trace stands for integration over momenta and conformal time. Differentiating with respect to the 
source Ja the generating functional Wa — —UnZ^, one can easily find the RG equations for the connected Green 
functions. The exact RG equation at any order in perturbation theory for the two-point correlation function is 
symbolically given in Fig. [5] where the dots indicated the fully renormalized vertex, the propagators are meant to 
be fully renormalized and the propagators with the square represents the RG kernel proportional to A-^{g^)^^ . 

RG methods are particularly suited to physical situations in which there is a separation between the scale where 
one is supposed to control the fundamental theory (in our case on subhorizon scales at the beginning of inflation) 
and the scale at which the measurements are actually made (on super horizon scales at the end of inflation). One 
can fix the initial conditions for a generic two-point correlation function Gpiji, T2) at the UV scale h-uv = {cipH). 
In this way, the correlator receives contributions only from those modes which are still inside the horizon, being 
ftp the scale factor at which the mode p exits the horizon, = (p/H). In this range of modes the correlator is 
supposed to be known and equal to the one given in Minkowski spacetime once the renormalization procedure is 
accounted for and all the parameetrs of the theory have been properly renormalized. Notice that, when computing 
the correlator for a given external momentum k, the UV cut off on the internal momenta will be chosen to be {atH), 
since one is interested in the effect of the internal superhorizon modes which are well out of the horizon when the 
mode k exits the horizon. 

The cut-off scale A is then decreased in such a way that more and more momenta corresponding to wavelengths 
outside the horizon are included. The cut off scale A is decreased till it reaches the smallest momentum which 
is excited during the de Sitter stage, that is up to A/^, — Auvio-i/cip) — H, where a,; is the scale factor at the 
beginning of inflation. THis corresponds to the comoving IR cut-off given by L^^. Starting from the fundamental 
scale, the RG flow describes the gradual inclusion of fluctuations at scales closer and closer to the one relevant to 
measurements. The key point is that the new fluctuations which are included at each intermeadiate step feel an 
effective theory which has been dressed by the fluctuations already included. 

The problem with this approach, however, is that the full set of RG equations for G"^ and has to be solved as 
both two-point functions are equally renormalized in the IR. Solving an approximate version of the RG equations 
corresponds to resum only a given class of diagrams, thus neglecting others which in most cases are not subleading at 
all. For instance, in the example developed in this Appendix of a cubic theory, neglecting the vertex renormalization 
and taking free propagators on the right-hand side of the RG equations amounts to resum that class of diagrams 
which at two-loop is represented by the first three graphs of Fig. 31 The inclusion of equally important diagrams 
in the IR amounts to solving the full set of RG equations with less and less approximations, making the approach 
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very complicated. 
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